function [outI, outQ, phase] = GSOP(Rxi,Rxq)

outI = Rxi/rms(Rxi);
outQ = Rxq - mean(Rxi.*Rxq)*Rxi/rms(Rxi)^2;
outQ = outQ/rms(outQ);

R1 = Rxi/rms(Rxi);
R2 = Rxq/rms(Rxq);

phase = asin(mean(R1.*R2));

tmp = -mean(R1.*R2);

aa = 1;


